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Abstract 

The Gaussian expansion method (GEM) is apphed to calculations of the nu- 
clear excitations in the random-phase approximation (RPA). We adopt the mass- 
independent basis-set that is successful in the mean-field calculations. The RPA 
results obtained by the GEM are compared with those obtained by several other 
available methods in Ca isotopes, by using a density-dependent contact interaction 
along with the Woods-Saxon single-particle states. It is confirmed that energies, 
transition strengths and widths of their distribution are described by the GEM 
with good precision, for the 1~, 2"'' and 3~ collective states. The GEM is then 
applied to the self-consistent RPA calculations with the finite-range Gogny DIS 
interaction. The spurious center-of-mass motion is well separated from the physical 
states in the E\ response, and the energy-weighted sum rules for the isoscalar tran- 
sitions are fulfilled reasonably well. Properties of low-energy transitions in ^'^Ca are 
investigated in some detail. 

PACS numbers: 21.60.Jz, 21.30.Fe, 21.10.Gv, 21.10.Pc 



Keywords: RPA calculation; Gaussian expansion method; finite-range interaction; giant 
resonance. 



1 Introduction 



As atomic nuclei far off the /3-stability exhibit exotic characters {e.g. nucleon halos [T] 
and new magic numbers [2]), excitation properties of unstable nuclei come under interest 
as well as their ground-state properties. The new features of unstable nuclei, such as 
the broad distribution of nucleons, the exotic shell structure and the coupling to the 
continuum, may give rise to new aspects of nuclear excitation. In theoretical studies of 
unstable nuclei, it is highly desired to handle those features efficiently. It is expected 
that ground-state properties of many nuclei can be investigated in the mean-field (MF) 
approximations, although correlation effects are not necessarily negligible. Among many 
efforts to develop new numerical methods to study structure of unstable nuclei, one of 
the authors (H.N.) proposed a method for MF calculations O IH [5] that is based on the 
Gaussian expansion method (GEM) [6]. The GEM is adaptable to the broad distribution 
of nucleons, even with finite-range effective interactions including semi-realistic nucleon- 
nucleon (NN) interactions [3, [8]. The results of the MF calculations suggest that the 
continuum effects are properly taken into account by the GEM [4]. Moreover, we have 
found that the basis parameters in this method are insensitive to nuclide, and thereby 
even a single set of bases can be applied to wide mass range of nuclei [5]. Because of these 
advantages, it is of interest to use the GEM for describing excitations of unstable nuclei. 

The random-phase approximation (RPA) provides us with a framework to treat one- 
particle-one-hole (Ip-l/i) excitations of nuclei in a consistent manner with the MF de- 
scription of the ground state. In this article we extensively apply the GEM to the RPA 
calculations. The strength functions of the spin-independent transitions are calculated 
in ^"'^^'^''Ca. We first use the single-particle (s.p.) states in the Woods-Saxon potential 
and adopt a density-dependent contact interaction for the residual interaction. For this 
schematic Hamiltonian, reliable results are obtained by the continuum RPA [9]. We test 
the new method by comparing its results with the continuum RPA results, particularly for 
excitations of unstable nuclei. There have been many RPA calculations using the harmonic 
oscillator (HO) s.p. basis functions. It is popular as well to implement RPA calculations 
by assuming the box boundary for the s.p. states. The new method is also compared with 
these methods. The present method is then applied to fully self-consistent Hartree-Fock 
(HF) plus RPA calculations, by employing the finite-range Gogny interaction. We shall 
check whether the spurious center-of-mass (cm.) excitation is well separated from the 
physical modes and how well the energy-weighted sum rules are satisfied. Properties of 
low-energy excitations in ^°Ca will also be discussed. 

2 Methods of calculations 

We shall consider the strength function, 

S{u;)=Y,6{u;-u;^) {a\O\0)\ (1) 



1 



where O stands for the transition operator, |0) the ground state, Ua the excitation energy 
of the state Below the particle-emission threshold we have discrete excited states. 
The sum over a in Eq. ([T]) is replaced by an energy integral for the continuum, if necessary. 
The response function is defined by 
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with r) representing an infinitesimal positive number. R{uj) is related to S{uj) through 
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By adding a finite imaginary number 27 to the energy u, we define 
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from which S{uj) is smeared out as 
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For the Ip-lh excitations S{uj) (or S^{uj)) is calculated from the transition amplitude 
(a|(9|0) in the RPA, or from R{uj) (or R^{uj)) in the linear response theory. It is well 
known that these two approaches are equivalent to each other [10]. We therefore do not 
distinguish them in this paper. 

We here discuss several methods of calculating S{ijj) or S^{uj) in the RPA. If we take a 
set of basis functions, the s.p. wave functions in the MF potential are given by superposing 
the bases. Solving the RPA equation that is represented by the s.p. states, we obtain 
the excited state |a) (to be more precise, forward and backward amplitudes along with 
cUq), from which ^(ti;) (or S^{uj)) is calculated. In this type of methods, quality of the 
results is governed by the s.p. basis functions. In the RPA calculations of nuclei near 
the /3-stability, it has been customary to employ the HO basis-set. However, it has been 
pointed out [3] that the HO set is impractical in handling the broad nucleon distribution 
and the coupling to the continuum, which could be important in unstable nuclei. In 
the present paper we propose a method using the GEM basis-set, which is comprised of 
multi-range Gaussian basis functions. Adaptable to the energy-dependent asymptotics of 
the s.p. wave functions including the oscillatory ones [1], the GEM basis-set is expected 
to describe the Ip-lh excitations appropriately, up to the continuum effects. 

It is also popular to obtain the s.p. wave functions by solving the MF equation in 
a discretized coordinate space. To keep the number of the s.p. states finite, one often 
confines the nucleons in a box or imposes a periodic boundary condition. In addition, 
truncation with respect to the s.p. energies is necessary in solving the RPA equation. 
While this method is feasible for contact NN interactions (or for local densities and 
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currents) such as the Skyrme interaction, it is not easy to be apphed when a finite-range 
interaction is adoptee^. 

In the hnear response theory R{uj) (or R^{uj)) is closely connected to the s.p. Green's 
function. The spectral representation, i.e. the expansion in terms of the s.p. states, 
may be employed for the s.p. Green's function [12]. It is also possible to take account 
of the exact form of the asymptotics of the s.p. Green's function at each uj, as was 
argued and carried out in Ref. [9]. This method is often called 'continuum RPA'. It is 
an additional advantage of the continuum RPA that we do not need truncation regarding 
the s.p. energies. The strength function S{uj) or S^{uj) is straightforwardly calculated 
from R{u!) or R^{uj), without constructing the excited state \a) explicitly. We note that 
the continuum RPA was extended to the quasiparticle RPA in Refs. [T3l [H]. Because 
it is implemented in a coordinate space, the continuum RPA is suitable for contact A^A^ 
interactions. 

There are other computational methods to implement the RPA calculations, although 
we do not treat them in this article. Recent developments include a method based on the 
mixed representation of the RPA equation [15], in which the unoccupied s.p. states are 
handled in a coordinate space [16]. The small amplitude limit of the time-dependent HF 
(TDHF) theory derives the RPA [10] as well. Methods to compute the linear responses 
via the TDHF have been explored [17]. One can carry out the TDHF calculations either 
by using s.p. bases or in a coordinate space. As a new approach based on the TDHF, the 
finite-amplitude method has been invented |18j . 

In all practical calculations in this paper, we assume that the ground states of the 
nuclei are well approximated by the spherical HF solutions. We shall compare results 
of several methods for ^°'^^'^°Ca in Sec. |3l In the following of this section we present 
computational details of the individual methods. We will focus on the method using the 
GEM in Sec. H 



2.1 Gaussian expansion method 

In the present paper we newly introduce a method of the RPA calculations using the 
Gaussian expansion method (GEM). In the GEM, both the ground and the excited states 
are represented by the s.p. basis functions having the following form: 

^.ijM = R.i,ir)[Y^'\v)x.]li^; 
Ruej{r) = Kij exp{-ur^) , (6) 

where Y^^\r) expresses the spherical harmonics and Xa the spin wave function. We drop 
the isospin index without confusion. The range parameter of the Gaussian basis function 
ly, which is also used as an index of the bases, can be complex (z/ = z/r -|- iu]) [6]. The 

^ It is mentioned that the Woods-Saxon basis functions, which arc obtained in a coordinate space, 
have been used in Ref. [ll] for a quasiparticle RPA calculation with the Gogny interaction. 
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constant Mutj is determined by 




(7) 



so as for {(puijml'^uejm.) to be unity. We take the set comprised of all the bases of 



with uq = (2.40 fm)~^ and b = 1.25 (6 real functions and 6 complex functions), irrespec- 
tively of With this basis-set we can describe ground states of many nuclei in a 
wide mass range from ^^O to ^osp]-, jg_ gj,g^ solve a MF equation in the s.p. space 
spanned by the above bases. We note that, since the GEM bases are not orthogonal to 
one another, the MF equation leads to a generalized eigenvalue problem. The MF solution 
determines the s.p. states, by which we represent the RPA equation. We then solve the 
RPA equation, following Sec. 8.4.4 of Ref. [lO]. 

2.2 Method of quasi-HO basis functions 

The HO basis-set has been employed for many MF and RPA calculations so far. The 
radial part of the spherical HO basis function has the form of a Gaussian multiplied 
by the associated Laguerre's polynomial. As mentioned in Ref. |3], we can produce a 
basis-set that is equivalent to the HO set by extending the basis functions of Eq. ([HD; 
R{r) oc r^'^'^P expl—ur"^) with an integer p, and u restricted to a single real value. While 
the individual basis functions are not the same as the HO functions, the set comprised 

of them is equivalent to the HO basis-set of {Lp^ ^\y/2i'r) r^exp(— i/r^);p = 0, 1, 2, ■ ■ ■ } 
for each where L^\x) is the associated Laguerre's polynomial, as revealed if we 

apply the Gram-Schmidt orthogonalization. An analogous basis-set is produced from the 
basis functions of Eq. ([6]), by taking u = b'"^^ (fc = 0, 1, 2, ■ ■ ■ ) with real and b close 
to unity. Indeed, if we take the b = 1 limit after orthogonalizing to the k = basis 

Ri,Qej{r) oc r^exp(— z/Qr^) = Lq^ ^\^/2l^r) expl—UQr^) , the k = 1 basis yields the HO 

basis of Lr^^ 

(■v/2z^r) r^exp(— i/q'"^)- It is easy to show that the basis proportional to 

L^^^ ^\y/2iyor) exp(— z/q?"^) is produced by successive orthogonalization, and thereby we 
obtain a set equivalent to the HO basis-set. Obviously 2k + i corresponds to the number 
of the oscillator quanta A'osc in the HO bases. In practice, we adopt a set of the basis 
functions of Eq. with b = 1.05 and Uq determined from cUqsc = 41.2A~^/^ MeV through 



uq = ^J2/ Mojosc- This set will be called 'quasi-HO' basis-set hereafter. The space spanned 
by the bases satisfying 2k +i < 11 is taken in the calculations in Sec. [31 With this basis-set 
we solve the MF and RPA equations in the same manner as in Subsec. 12. 1[ 




(A; = 0,1,--- 

(A; = 0,1, 2) 
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2.3 Coordinate-space method with box boundary 

If the MF is local, it is not difficult to solve the MF equation in a coordinate space. 
Because we assume the spherical symmetry, it is sufficient to consider one-dimensional 
coordinate space of the radial variable. This radial coordinate is discretized by a mesh 
whose size is denoted by h. We impose a box boundary condition that the s.p. wave 
functions should vanish at rmax- The MF equation is solved in the coordinate space thus 
determined. In Sec. |3]we adopt h = 0.2 fm and rmax = 20 fm, conffiming the convergence 
in the MF calculations. 

The RPA Hamiltonian is represented by the s.p. states obtained in the MF calculation, 
by cutting off the unperturbed excitation energy at Ascut- In practical calculations in 
Sec. 131 we do not solve the eigenvalue problem defined by the RPA equation, but we instead 
compute R^{uj) by solving the Bethe-Salpeter equation in the spectral representation [12] 
(using the s.p. states and Zi^cut)- The strength function S.^{ijj) is then obtained via Eq. ([3]) 
(or ([5])). For certain transition modes the RPA results are slow to converge for increasing 
^£cut [in]- We have con&med that S^{(jj) at < 30 MeV has no visible difference between 
^^cut = 150 MeV and 200 MeV for all the modes under consideration. 

2.4 Continuum RPA 

In the continuum RPA [9] , the strength function is calculated from the response function. 
However, the response function is constructed from the s.p. orbits of holes and the 
s.p. Green's function for the excited particles, unlike the spectral representation. The 
asymptotic form of the s.p. Green's function at large r, which depends on cu, is known for 
each partial wave. For each cu, the inner part of the s.p. Green's function is continuated 
to the proper asymptotic form at a sufficiently large r. Hence the calculation is free both 
from the boundary at Tmax and from the energy cut-off like Zi^cut- For a relatively simple 
interaction as used in Sec. [3l we can obtain 'exact' transition strength functions (which 
are exact within the RPA), as long as the convergence with respect to h is reached. In the 
continuum RPA calculations in Sec. [3l we follow the method shown in Ref. [13], except 
that we do not need integration in the complex energy plane because we do not have 
pair correlations in the ground state. The mesh size h is taken to be 0.2 fm as in the 
box-boundary calculation, for which the convergence has been conffimed. 

Whereas the continuum RPA calculations are implemented with specifying uo, it is 
not easy to pinpoint the excitation energies of the discrete states. To avoid missing the 
discrete levels, we calculate S^{uj) rather than S{uo), taking 7 = 0.2 MeV. The excitation 
energy and the transition strength of each discrete state is extracted by ffiting S^{uj) 
around the peak to a Lorentzian. Correspondingly, we compute S^{uj) with 7 = 0.2 MeV 
in the other methods, by using the imaginary energy + ^7 in the method of Subsec. 12.31 
and by smearing out S{uj) in those of Subsecs. 12. II and 
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3 Comparison of methods for contact force 



In this section numerical results of the four methods explained in the preceding section 
I2.1H2.4I are compared. For this purpose, we take the nuclear MF of the Woods-Saxon 
potential to which the Coulomb potential is added, as adopted in Ref. [9]: 

U{r) = (l- 0.67^^r,) (uo f{r) + U,,e-s l^f{r?j + ^(1 - r,) Uc{r) , (9) 

where /(r) = 1/[1 + exp((r - R)/a)], Uo = -58MeV, Uis = SOMeVfm^, R = 1.20(A - 
ly/^im, a = 0.65 fm, = +1 (—1) for a neutron (a proton), and Uc{r) is the Coulomb 
potential produced by the uniform charge distribution of [Z — 1) protons in the sphere of 
radius R. For the residual interaction, we assume the density-dependent contact force as 

^'res = f\to{l + XoP^) 5(ri - Fa) + ^^3(1 + XsP^) p(ri) 5(ri - ra)! , (10) 

L 5 J 

with to = -llOOMeVfm^ xo = 0.5, = 16000 MeVfm*^ and X3 = 1 0. p(r) denotes the 
nucleon density, and Po- the spin exchange operator. It will be stated below how to fix the 
overall factor /. The spin densities arising from gnored for the sake of simplicity. 

We shall consider transitions carried by the one-body operator, 

with A = 1, 2, 3, where i denotes the index of nucleons, for the ^°Ca, ^^Ca and ^°Ca nuclei. 
The (A = 1, r = 0) mode corresponds to the spurious cm. motion. Subtracting its 
contribution, we obtain a modified operator for the (A = 1, r = 1) mode, 

which is proportional to the El operator with the cm. correction. The expression i G n 
{i G p) indicates that the sum includes the i-th nucleon if it is a neutron (proton). The 
strength functions S^-^''^^ and S'^^'^-' are defined by substituting O^'^''^^ in Eqs. (fT][5l) . For 
the (A = l,r = 1) mode we can define ^(^=1.^=1) and 5^^=^'^=^) from (9(^=1-^=1). The 
renormalization parameter / in Eq. ffTOl) is determined so that the spurious cm. state 
should have zero excitation energy [5| , for individual cases. We tabulate the / values thus 
determined in Table [TJ Although there is no consistency between the MF and the residual 
interaction, this set of the MF potential and the interaction is suitable for assessing the 
numerical methods. 

In Figs. [TH5]we display the strength functions 5'^'^'^-' (cj) . For the calculations assuming 
the box boundary shown in Subsec. 12. 3[ we take Ascut = 50 MeV, whose results are not 
qualitatively different from those of ZiScut = 200 MeV. 

If we estimate from the s.p. energies, the particle threshold lies at u = 8.7 MeV for 
^°Ca and ^^Ca, while at 3.4 MeV for ^''Ca. In all the three nuclei we find A = 3 peaks 
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Table 1: Adopted values of /, the renormalization factor of the residual interaction. 
The label 'Cont.' ('Box') indicates the case of the continuum RPA (the box-boundary 
method). 



nuclide 


Cont. 


Box 


quasi-HO 


GEM 


40Ca 


0.691 


0.785 


0.715 


0.700 


^^Ca 


0.711 


0.808 


0.750 


0.720 


60Ca 


0.758 


0.842 


0.796 


0.764 



at low uj that correspond to discrete states. For ^^Ca we also find low- lying discrete 2"*" 
states because of the shell structure. Although the 8.86 MeV state in the A = 3 excitation 
of ^^Ca is located just above the neutron threshold, it behaves like a discrete state. For 
these discrete states we list the excitation energies Ua and the transition strengths 



(13) 



in Table O The low-lying 3~ state is highly collective in any of the three nuclei. For this 
state the convergence for Zi^cut is quite slow in the calculation under the box boundary 
condition, and we show the results with Ascnt = 200 MeV in Table El Note that the 
5'^'^''^)(ci;) graphs in Figs. [Hand [5] do not much change even if we take Ascut = 200 MeV. 



7 




5 10 15 20 25 



20 n 




5 10 15 20 25 30 

(0 [MeV] 



Figure 1: S'^'^"^'^"^^ (w) in ^°'^^'^°Ca, by the RPA calculations using the mean field of Eq. ([9]) 
and the residual interaction of Eq. ffTOl) . Results of the continuum RPA, the box-boundary 
method, the method employing the quasi-HO bases and the GEM are represented by the 
black solid, the black dashed, the green solid and the red solid lines, respectively. We take 
7 = 0.2 MeV for all the calculations. 
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Figure 2: S!^^=^'^=°\cu) in ^^'^^'^^Ca. See Fig. □ for conventions. 
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Figure 3: S!,^=^'^=^\cu) in ^"'^S'^^Ca. See Fig. □ for conventions. 
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Figure 4: S!,^=^'^=°\cu) in ^"'^S'^^Ca. See Fig. □ for conventions. 
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Figure 5: S^^~^''^~^\oi!) in ^"'^'^'^^Ca. See Fig. [T] for conventions. 
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Table 2: Excitation energy uOa (MeV) and transition 
strength B^''^^ (fm^-^) for each discrete state, obtained by 
the RPA calculations employing the mean field of Eq. ([9]) 
and the residual interaction of Eq. (fTOl) . Results of the 
continuum RPA, the box-boundary method, the method 









using the quasi- 


•HO bases and the GEM are 


compared. 








nuclide 








Cont. 




Box 


quasi-HO 




GEM 






A 


r 




5(A,r) 

a 




a 






UJa 


5(A,r) 
a 


40Ca 


3 





2.96 


8.49 X 10"^ 


3.04 


8.20 X 10^ 


3.28 


7.51 X 10^ 


3.00 


8.32 X 


10^ 






1 




3.29 X 10^ 




3.25 X 102 




2.97 X 102 




3.20 X 


102 


^^Ca 


2 





2.45 


7.88 X 102 


2.45 


7.87 X 102 


2.49 


7.78 X 102 


2.45 


7.89 X 


102 






1 




2.03 X 102 




2.02 X 102 




1.99 X 102 




2.03 X 


102 




2 





5.92 


7.38 X 10^ 


5.92 


7.38 X 10^ 


5.93 


6.68 X 10^ 


5.92 


6.85 X 


101 






1 




1.57 X 10^ 




1.56 X 10^ 




1.42 X 10^ 




1.43 X 


101 




3 





5.18 


5.82 X 10^ 


5.20 


5.77 X 10^ 


5.27 


5.66 X 10^ 


5.19 


5.77 X 


10^ 






1 




1.33 X 10^ 




1.36 X 10^ 




1.41 X 10^ 




1.31 X 


10^ 




3 





8.86 


2.22 X 10^ 


8.86 


2.24 X 10^ 


8.89 


2.17 X 10^ 


8.86 


2.17 X 


10^ 






1 




5.35 X 10^ 




5.48 X 10^ 




4.68 X 10^ 




4.81 X 


10^ 


6°Ca 


3 





1.96 


4.07 X 10^ 


1.98 


4.02 X 10^ 


2.09 


3.91 X 10^ 


1.97 


4.06 X 


10^ 






1 




1.15 X 10^ 




1.15 X 10^ 




1.17 X 10^ 




1.15 X 


10^ 



By comparing with the results of the continuum RPA, we find that the maximum 
deviation of uo^ is 0.08, 0.32 and 0.04 MeV in the box-boundary, the quasi-HO and the 
GEM results, respectively. The relative error of B^''^^ is 3% in the box-boundary, 13% in 
the quasi-HO and 10% in the GEM results, at maximum. Thus, as long as the discrete 
states are concerned, the GEM provides precision of ~ 0.05 MeV for Ua and ~ 10% for 
Ba- The precision of the GEM for is slightly better than that of the box-boundary 
calculations, and substantially better than that of the calculations using the quasi-HO 
bases. Although the precision might look worse than the box-boundary method for B^, 
the strong transitions are described with good precision. It is noticed that the quasi-HO 
basis-set gives rise to sizable errors for the lowest-lying 3~ states. 

Above the particle-emission threshold, the continuum RPA gives smooth S{uj) {i.e. 
the 7—^0 limit of S^{uj)), whereas in the other methods S{uj) is still represented by 
a sum of the delta functions. In Figs. [THS] the strength functions are smoothed by the 
width parameter 7 (= 0.2 MeV). It is a clear advantage of the continuum RPA that the 
smooth behavior of 5'(a;) in the continuum is automatically taken into account, and could 
be important in investigating certain aspects of transitions, e.g. the soft dipole responses 
near the threshold [131 120] • However, at higher u the strength functions should further 
be smeared due to coupling to the 2p-2h degrees of freedom, which is not incorporated 
in the usual RPA. Therefore gross properties of the strength distribution will be more 
important than its fine structure, when we view the RPA results. In Table [3], we tabulate 
the transition strength, the average energy and its standard deviation in an energy domain 



The domain D is chosen so as to cover the broad resonance- like structure of S^'^''^\uj). 
For the (A = 1, r = 1) mode, we use S!^^^^''^^^\uj) instead of S!^^^^''^^^\uj) as before. 




(14) 
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Table 3: Comparison of B^^'^' (fm^^), oo)^'^> (MeV) and 
(T^'^"* (MeV) in a certain energy range D, which is spec- 
ified by uj'B"' and cug^^ (MeV). 



nuclide 














Cont. 








Box 








quasi 


-HO 








GEM 






A 


r 




D 








^(A,r) 
D 






" D 




r) 




" D 




-,(A,r) 
D 


^(A,r) 
D 




1 


1 


10.0 


- 25.0 


7 78 


X 


10" 


17.35 


2.38 


7 /I/I 
( .44 


X 


iU 


17.85 


2.40 


7 79 


X 


lo'J 


17.54 


2.38 


7 S9 
( .0/ 


X 


iU 


17.43 


2.43 




2 





12.0 


- 22.0 


i.DD 


X 


10^ 


17.12 


1.10 


1 fi7 
i.D / 


X 


1 n3 

iU 


16.88 


1.14 


i.Do 


X 


10^ 


17.23 


1.12 


i.DD 


X 


1 n3 

iU 


17.09 


1.10 






1 


12.0 


- 30.0 


D.Di 


X 


102 


24.49 


3.98 


o.yo 


X 


iU 


24.29 


3.99 


D.4o 


X 


102 


24.65 


3.72 


D.oo 


X 


in2 


24.73 


4.06 


48Ca 


1 


1 


10.0 


- 25.0 


9.33 


X 


10° 


17.06 


2.45 


8.92 


X 


10° 


17.61 


2.47 


9.29 


X 


IQO 


17.37 


2.56 


9.40 


X 


10° 


17.15 


2.49 




2 





12.0 


- 22.0 


2.29 


X 


10^ 


16.32 


1.16 


2.28 


X 


10^ 


16.13 


1.13 


2.26 


X 


10^ 


16.30 


1.07 


2.29 


X 


10^ 


16.30 


1.19 






1 


12.0 


- 30.0 


1.08 


X 


10^ 


23.36 


4.67 


9.90 


X 


102 


23.09 


4.80 


1.04 


X 


10^ 


23.40 


4.70 


1.09 


X 


10^ 


23.43 


4.70 




3 


1 


10.0 


- 22.0 


1.71 


X 


10^ 


15.08 


3.22 


1.62 


X 


10^ 


15.30 


3.25 


1.55 


X 


10^ 


14.65 


3.20 


1.70 


X 


10^ 


15.11 


3.27 


60Ca 


1 


1 


2.0 


- 10.0 


2.11 


X 


10° 


6.75 


1.54 


1.90 


X 


10° 


6.71 


1.55 


1.84 


X 


10° 


7.17 


1.07 


2.09 


X 


10° 


6.81 


1.49 








10.0 


- 25.0 


1.04 


X 


10^ 


16.24 


2.86 


1.01 


X 


IQi 


16.74 


2.95 


1.04 


X 


IQi 


16.54 


2.89 


1.04 


X 


10^ 


16.31 


2.84 


h- ' 

OX 


2 





10.0 


- 20.0 


4.09 


X 


10^ 


14.48 


1.52 


4.09 


X 


10^ 


14.38 


1.53 


4.18 


X 


10^ 


14.24 


1.56 


4.09 


X 


10^ 


14.42 


1.54 






1 


3.0 


- 12.8 


8.85 


X 


10^ 


9.06 


2.32 


8.72 


X 


102 


8.98 


2.25 


8.78 


X 


102 


10.80 


1.43 


8.92 


X 


102 


8.98 


2.16 








12.8 


- 20.0 


1.28 


X 


10^ 


15.16 


1.54 


1.31 


X 


10^ 


15.01 


1.47 


1.11 


X 


10^ 


15.25 


1.44 


1.30 


X 


10^ 


15.30 


1.43 




3 





3.5 


- 15.0 


1.30 


X 


10^ 


7.97 


2.72 


1.31 


X 


10^ 


7.91 


2.73 


8.66 


X 


10^ 


7.89 


2.66 


1.34 


X 


10^ 


7.94 


2.67 








15.0 


- 30.0 


8.49 


X 


10^ 


23.78 


4.00 


8.35 


X 


10^ 


23.74 


3.91 


9.84 


X 


10^ 


22.49 


3.85 


8.91 


X 


10^ 


24.09 


3.95 




3 


1 


3.5 


- 15.0 


1.31 


X 


10^ 


8.47 


2.86 


1.28 


X 


10^ 


8.48 


2.89 


8.15 


X 


10^ 


9.08 


2.98 


1.32 


X 


10^ 


8.43 


2.80 








15.0 


- 30.0 


6.13 


X 


10^ 


22.81 


4.15 


6.13 


X 


10^ 


22.88 


4.12 


7.73 


X 


10^ 


21.73 


4.11 


6.40 


X 


10^ 


23.07 


4.19 



In Table [3] we find that the errors in the GEM results are 0.3 MeV for uJd''^ , 5% for 
B^]^'^^ and 0.2 MeV for (T^'^"* at maximum, which are compared to 0.5 MeV (1.1 MeV), 
10% (30%) and 0.1 MeV (0.9 MeV) in the box-boundary (the quasi-HO) resuhs. We 
thus confirm that the GEM reproduces distribution of the transition strengths with good 
precision up to from stable to drip- line nuclei. 

In Figs.lIHSl we view considerable transition strengths at low u in ^°Ca, which arise due 
to the loosely bound neutrons. Similar results have been reported in the continuum RPA 
calculations with the Skyrme interaction [211 [22l [23] . Whereas such low-lying transition 
strengths in drip-line nuclei will be investigated further in a forthcoming paper, we here 
argue dependence of the results on the computational methods. In comparison to the 
results obtained with the other methods, notable deviation is found in the quasi-HO 
results for the strengths at < 10 MeV. This is because the quasi-HO basis functions 
are not suitable for describing loosely bound nucleons which may have broad spatial 
distribution. To clarify this point, we consider the transition density defined by 

rVgUr;«) = {a\Y.6{r~r,)r^Y^^\rm {T,=p,n). (15) 

In Fig. [HI the neutron transition densities are depicted for the low-energy A = 3 modes 
of ^°Ca. We present the transition densities at two prominent peaks in 5 < u; < 7.5 MeV 
obtained in the GEM and the quasi-HO calculations. For comparison, the densities in 
the continuum RPA are also displayed, whose renormalization factors {C in Eq. (17) of 
Ref. [20]) are determined so that the densities should be comparable to their counterparts 
in the GEM and the quasi-HO results. 

It is found that the transition densities depend primarily on the s.p. bases, i.e. the 
GEM or the quasi-HO, rather than on the excitation energies of the states. These low- 
energy transitions are dominated by excitation of a neutron from the p/-shell to the 
continuum. The rapid decrease of r'^p[^~^\r;a) at r ^ 9fm in the quasi-HO results is 
attributed to limitation of the basis functions mentioned above. It does not seem easy 
to describe such transitions properly, even if we increase the number of the quasi-HO 
basis functions. On the contrary, the GEM gives broad distribution of p[^~^\r; a) , in 
fair agreement with the corresponding transition densities in the continuum RPA. This 
implies that by the GEM we can take into account the effects of coupling to the continuum 
sufficiently. We here comment that r'^ Ptr^'^\f", «) obtained with the box-boundary method 
is also close to that of the continuum RPA in this energy region. 

We have thus established that the method using the A-independent GEM basis func- 
tions of Eq. ([8]) is capable of describing the excitation energies, the transition strengths 
and the widths of their distribution with good precision, including nuclei in vicinity of 
the drip line, within the RPA framework. This method gives precision comparable to the 
method assuming the box boundary, and substantially better than the method using the 
quasi-HO basis functions for nuclei near the neutron drip line. The present method based 
on the GEM cannot reproduce the smooth energy dependence of the transition strengths 
in the continuum region. This is a defect if compared to the continuum RPA method. On 
the other hand, the GEM has an advantage in adaptability to finite-range interactions. 
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Figure 6: r'^Pir,n {r^a) for low-energy states a in ®°Ca. The red solid (dashed) line is 
obtained by the GEM at u;„ = 5.02 MeV (6.68 MeV), while the green solid (dashed) 
line by the quasi-HO basis functions at oj^ = 5.11 MeV (7.37 MeV). All these states carry 
relatively strong transition strengths of B\^~^''^^^^ > 5 x 10^ fm^. The black solid (dashed) 
line is the continuum RPA resuh at u; = 5.02 MeV (6.92 MeV). 

In the next section we shall apply the present method to self-consistent RPA calculations 
with a finite-range interaction. 

4 Application of GEM to self-consistent RPA 

The numerical method based on the GEM is applied to self-consistent RPA calculations. 
We consider the effective nuclear Hamiltonian comprised of the kinetic energy and the 
effective NN interaction, 

Hr, = K + V^- K = Y.^, V5v = E%- (16) 

i ^^^^ i<j 

Here i and j are indices of the constituent nucleons. The full Hamiltonian is given by 
H = if TV + Vc~ Hera., where Vc stands for the Coulomb interaction among protons and 
Hem. is the cm. Hamiltonian. We use the A-independent GEM basis functions of Eq. ([8]) 
for all the calculations in this section. Once we compute and store the two-body matrix 
elements of Vat, we can use them both for the MP and the RPA calculations (except the 
density-dependent part of Vjv). The exchange term of Vc is exactly treated [5]. It should 
be commented that Hem. is highly non-local, and therefore is not easy to be handled in 
the coordinate representation. In contrast, we can easily take into account both the one- 
and two-body terms of i^c.m. in the GEM. 

In the calculations shown below, we adopt the DIS parameter-set [24] of the Gogny 
interaction [25] for Vn, in which the central channel has finite ranges, and take ^°'^^'^°Ca 
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Table 4: lu"^ values (MeV^) in the self-consistent RPA calculation by the GEM with the 
DIS interaction. 



nuclide 




40Ca 


-5.80 X 10^6 


48Ca 


-8.61 X 10-6 


60Ca 


-2.67 X 10-6 



as examples. Note that ^^Ca is located near the neutron drip line in the prediction with 
DIS [8]. We carry out the RPA calculations on top of the spherical HF solutions for these 
nuclei. 

In Ref. [26], the DIS interaction was applied to the quasiparticle RPA calculations 
for the Si and Mg nuclei near the /3-stability line, by employing the HO basis functions. 
The present method based on the GEM will be useful for extending such approaches to 
drip-line nuclei. 



4.1 Spurious cm. motion 

It is proved that the spurious states, which are the Nambu-Goldstone modes emerging 
from the spontaneous symmetry breaking, have zero excitation energy and are separated 
from the other states in the self-consistent RPA [10]. A typical example is the spurious cm. 
motion. However, we do not necessarily have a zero-energy state in practical calculations 
even if the Hamiltonian keeps the translational invariance, because of finite size of the s.p. 
space. Conversely, the energy of the spurious cm. state provides a measure of accuracy 
of the numerical calculation. 

We here express the lowest 1- state in the solution of the RPA equation, which cor- 
responds to the spurious cm. motion, by a = s. In Table Hj Ug values obtained by the 
GEM basis-set of Eq. ([8]) are shown for ^"'^^'^''Ca. By the GEM the spurious state has zero 
energy to good precision. Moreover, if the spurious state is well separated, the transition 
strength 5^^=^'^=°) should vanish for a ^ s. We confirm B^^^''^^'^^ < 10-'^ fm^ for all 
a{y^ s), in any of the three nuclei. 

4.2 Energy- weighted sum rules 

The RPA preserves the energy- weighted sum rule (EWSR), 

Si = Y.uJ{a\O\0)f = l{0\[O\ [H, O]]\0) , (17) 

if the expectation value of the double commutator is evaluated for the MP ground 
state [lO]. In correspondence to the specific mode O*^^'"^-', we define S^^'^''. The energy- 
weighted sum S^'**'^'' is associated with the quantities in Eq. (fT4l) by S^'^'^'* = ci)^'^^ ■ -B/)''^^ 
for D covering the whole range of lo {i.e. — u^'^^ = oo). In evaluating 
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the double commutator, we first take H' = + Vc, ignoring Hem., instead of tlie full 
Hamiltonian H. Most effective interactions (or energy density functionals) including the 
Skyrme, the Gogny and the M3Y-type ones [71 [8] lead to the continuity equation for the 
isoscalar density and current (as long as we use H'), from which the following EWSR for 
the isoscalar transition is derived [27], 



^/(A,r=0) 



Lo|[o(^.-=o)tj/7',0(A,-=o)]]|o) 



A(2A+1)2 1 



2M 



(0|E 



2A-2 



|0) 



where we restrict ourselves to A > 2. Though not written explicitly, the z-components of 
are summed up. However, ifc.m. modifies the EWSR as 



El' 



(A,r=0) 



Lo|[C)(^.-=o)tj/^^C'{A,-=o)]]|o) 
A(2A+1)^ 1 r^o|Erf-|0) 



47r 2M 

47r 1 / _ 

-2a3TI<°I(E 



,A-ly{A-l) 



A(2A + 1)2 1 



2M 



1 



A 



(oiE 



2A-2 



|0) 



2A 



(19) 



For the A = 2 case, Eq. (|T9|) can be rewritten as 



(A=2,T=0) 



50 1 



47r 2M 



(OlE(r.-RflO), 



(20) 



where R denotes the cm. position. This correction due to Hem. is harmonious with 
the cm. correction to the rms matter radius [7]. The EWSR of Eq. f|T9|) can be used 
as another tool to check appropriateness of the s.p. space adopted in the numerical 
calculation. 

For the (A = l,r = 1) mode, the Thomas-Reich-Kuhn (TRK) sum rule is derived if 
Vn is local. 



Strk = lmd^'=''^='^\ [K, 0(^=1'^=^)]] |0) = 



9 1 AZN 



4:71 2M A 



(21) 



However, S^'^ ^'^ shifts from Strk in practice, owing to the non-locality in the charge 
exchange term of Vat. The enhancement factor k is defined by 



El- 



(A=l,r=l) 



:i + /«)s 



TRK 



(22) 
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It is noted that -ffc.m. does not influence the energy-weighted sum of the (A = 1, r = 1) 
mode {i.e. [O^^^^'^^^^'^ , [Hem., O^'^^^''^^^^]] = 0), as long as it is fully taken into account. 
The enhancement factor k, can be expressed by an expectation value of the charge exchange 
term of Vn for the ground state. In the case of the Skyrme interaction, an explicit 
expression of k in terms of the neutron and proton densities has been given in Ref . [28] . 
It is complicated to evaluate k in an individual nucleus if V^- is taken to be a finite-range 
interaction. However, k in the symmetric nuclear matter is related to the Landau-Migdal 
parameter /{ [7] by 

where M* denotes the effective /c-mass. Equation fl23l) gives an approximate value of k 
for finite nuclei. 

In many calculations using the Skyrme interaction, only the one-body term of i^c.m. 
is taken into account. We point out that this prescription influences the above sum rules. 
For the (A > 2, r = 0) modes, the energy-weighted sum becomes multiplied by 

(1 — 1/A), lacking the second term in the last expression of Eq. (fT9|) . Moreover, if we 
subtract only the one-body term of Hem., the TRK part of the energy- weighted sum for 
the (A = l,r = 1) mode {i.e. Eq. ( j2Tl) ) is fictitiously reduced by (1 — 1/A). 

As an indicator for the EWSR, we define the ratio, 

7^{^'-) = ^u;^ (a|O(^'")|0) ' /^(0|[O(^'")^[i^,O(^'")]]|0), (24) 

a I ^ 

which should be unity if the s.p. basis-set is complete. In Table [5], we show 7^(-^''^=°) 
(A = 2, 3) in ^°'^^'^°Ca. It is confirmed that the present method satisfies the EWSRs of 
Eq. (|T9l) within a few percent precision. For the (A = 1, r = 1) mode, we use the ratio to 
the TRK sum rule, 

^(A=i,.=i) ^ 2 [K,a(^=i'"=^)]]|0) , (25) 

a / 2 

and compare it with the rhs of Eq. (123|) . which is denoted by '^xoo^"^'*- The values 
of :^SJ=^'^=^^(= 1 + k) are presented in Table O The ^'^ value does not depend 
strongly on nuclides. The nuclear matter value of 'f^K~^^~^^ gives a good first approxima- 
tion, and slight reduction from "^xoo'^"^^ accounted for mainly by the difference 
between Pp(r) and p„(r) [2B]. It is noted that the experimental value has been reported 
as ■^^"^'^"^^ = 1.76 ± 0.10 [29j, being almost independent of the mass number. 



4.3 Strength functions 

We next show the strength functions S[^''^^{uS) obtained by the self-consistent RPA cal- 
culations with DIS, in Figs. [Tl-fTTl We here adopt 7 = 0.4 MeV. The strength functions 
in the RPA are compared with the unperturbed strength functions, in which the residual 
interaction is ignored. 
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Table 5: 7^(-'*'^-o) in the RPA calculation by the GEM, in which the DIS interaction is 
employed. 

^iuclid^^~^T^^^^ 



40Ca 2 1.005 

3 1.031 

^^Ca 2 1.006 

3 1.033 

^OQa 2 1.003 

3 1.010 



Table 6: n^i~^'^~^^ in the RPA calculation by the GEM, in comparison with the nuclear 
matter value UkTc^'^^^^ = (M/M*) (1 + f[/3). The DIS interaction is used. 



nuclide 


^(\=l,T=l) ^(A=l,r=l) 


^oCa 


1.587 


^^Ca 


1.589 1.660 


60Ca 


1.584 



Irrespectively of (A, r), the strength functions S!:^'^''^\lu) with the DIS force are qualita- 
tively similar to those obtained with the schematic interaction shown in Sec. O Significant 
difference between the RPA strength function and the unperturbed one suggests collec- 
tivity of the transition. In this regard we view collectivity of the isovector giant dipole 
resonance (IV-GDR) induced by a repulsive part of the residual interaction in Fig. [71 and 
of the isoscalar giant quadrupole and octupole resonances (IS-GQR and IS-GOR) induced 
by an attractive interaction in Figs. M and [TUl Collectivity of these giant resonances is 
further confirmed from the forward and backward amplitudes of the state that forms the 
highest peak in the strength function; a number of the unperturbed excitations are mixed 
in the RPA states. 

In ^°Ca, transition strengths emerge and form a peak in S'^'^''^^ at relatively low energy 
(ti; < 10 MeV), though the transitions are not so strong as in the giant resonances. It is of 
interest whether such low-energy strengths have collective nature or not. By comparing 
the RPA strength functions to the unperturbed ones, it is suggested that the low-energy 
transitions specific to ^°Ca hardly have strong collectivity. This consequence is consistent 
with the argument in Ref. [30], and will further be investigated in a forthcoming paper. 
For the (A = 1, r = 1) mode, the low-energy strengths may be compared to the so-called 
pygmy dipole resonance (PDR) that has been observed in several Z < N nuclei such 
as ^^Zt [31] and ^osp^ fj^j^g IV-GDR in Z ~ nuclei is interpreted as a surface 

oscillation to which the protons and the neutrons contribute with opposite phases. In the 
neutron-rich nuclei the El transition could have two collective components; an oscillation 
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Figure 7: S\^=^'^=^\uj) in ^o.48,60q^_ ^^j^g ^ ^ 0.4 MeV. The solid lines are obtained 
by the self- consistent RPA calculations with the DIS interaction, while the dashed lines 
represent the unperturbed strength functions. The particle threshold energy in the HF 
approximation is indicated by the inverted triangle for each nucleus. 
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Figure 8: S'^'^ ^'"^ ^K^) in ^°'^^'^°Ca calculated with the DIS interaction. Conventions are 
the same as in Fig. [71 
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Figure 9: S'^'^ ^'"^ ^K^) in ^°'^^'^°Ca calculated with the DIS interaction. Conventions are 
the same as in Fig. [71 
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Figure 10: 5'^'*' ^'"^ "H^) ^°'^^'^°Ca calculated with the DIS interaction. Conventions 
are the same as in Fig. [71 
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Figure 11: 5'^'*' ^'"^ ^H^) ^°'^^'^°Ca calculated with the DIS interaction. Conventions 
are the same as in Fig. [71 
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between protons and neutrons, and an oscillation between a core and excess neutrons. 
The latter is a possible interpretation of the PDR. In Fig. [12] the neutron and proton 
transition densities in ^''Ca are depicted for the uja = 6.9 MeV and 16.5 MeV states, both 
of which form the peaks in Fig. [71 We here define the transition density by 



r;a) 



n] 



' 22' 

— i^a\Y^b[T-r,)r,Y^^\im (for r, 

— (a| ^5(r - r,) r,y«(f,)|0) (for = p) 



(26) 



which is analogous to Eq. f[T^ but is subject to the modification of Eq. f[T^ . Because the 
spurious cm. motion is well separated, we have 



where the isoscalar transition density is defined by 



2 (A=l)/ N 



A 



2NZ 



N ■ r^plr,n^ (r; a) + Z ■ r^pl'^^p '[r;a : 



„2~(A=1), 



(27) 



:28) 



Figure [T2l shows that, to the first approximation, the higher-lying strength corresponds to 
an out-of-phase oscillation between protons and neutrons as in the usual IV-GDR, while 
the lower-lying strength looks like an oscillation between a core and outer neutrons, as 

2~(A=1), 



pointed out in Ref. [33]. It should be noted that the peaks in r plj.~^ '{r; a) are considerably 
displaced between = n and p for the transition to the Ua = 16.5 MeV state. This 
indicates that, whereas r'^ p[^Ys\''^'i vanishingly small in the Z ^ N nuclei such as 
^°Ca, it does not vanish in ^°Ca even for the strong transition, though satisfying Eq. ([27|) . 
Moreover, behavior of r'^pt^^\r:, a) at r 4fm suggests that this state contains a weak 
admixture of the oscillation between a core and outer neutrons. 



5 Summary 

We have developed a method of implementing RPA calculations, which is based on the 
Gaussian expansion method (GEM). Owing to the advantages of the GEM which have 
been established in the MF calculations, it is naturally expected that we can efficiently 
compute excitation of nuclei including coupling to the continuum. The parameters of the 
s.p. basis functions are insensitive to nuclide, and even calculations with a single set of 
bases may cover wide range of the mass table. 

The method has first been tested in '^"'^^•^^Ca with a density-dependent contact inter- 
action on top of the Woods-Saxon single-particle (s.p.) states, by comparing its results 
with the results obtained with the continuum RPA method. We have confirmed that 
the present method using the GEM basis functions of Eq. ([8]) describes the energies, the 
transition strengths and the widths of their distribution with good precision for the 1~, 
2+ and 3~ collective states, including drip-line nuclei, although we cannot reproduce the 
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Figure 12: r'^p\^ .^J{r]a) in ^'^Ca for the Ua = 6.9 and 16.5 MeV states, obtained by the 
RPA calculation with the DIS interaction. The blue (red) line is for = n (p). 
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smooth energy dependence of the transition strengths in the continuum. The new method 
has been compared also with the box-boundary method and with the method using bases 
analogous to the harmonic-oscillator (HO) ones. The present method attains precision 
comparable to the box-boundary method of rmax = 20 fm. The basis-set similar to the 
HO set of iVosc < 11 does not give good precision, particularly for low-energy transitions 
in drip-line nuclei, because it is unable to reproduce transition densities at r > 10 fm. It 
may not be easy to describe such transitions appropriately even if we take a larger number 
of the HO basis functions. 

Another advantage of the present method is its tractability of finite-range interactions. 
This point has been demonstrated by the application to the self-consistent HF plus RPA 
calculations in ^'''^^'^"Ca with the Gogny DIS interaction. It has been confirmed that 
the zero excitation energy of the spurious state and the energy-weighted sum rules for 
the isoscalar transitions are fulfilled to good precision. By comparing the RPA strength 
functions to the unperturbed ones, collectivity of the transitions has been argued. For 
^°Ca, characters of the low-energy dipole strengths have been investigated as well as those 
of the giant dipole resonance, via the transition densities. 

Further application of the present method to excitations of nuclei with finite-range 
interactions, including the semi- realistic interactions [3 |8], is under progress. We here 
mention a study of the Ml transition in ^osp]^ in which significant effects of the 
tensor force have been confirmed. Future plans include extension of the present method 
to the quasiparticle RPA, which is quite promising since the GEM has been successfully 
applied to the Hartree-Fock-Bogolyubov calculations [1]. 
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